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Abstract 

It is shown that non-adiabatic cosmological perturbations cannot appear 
during the period of reheating following inflation with a single scalar inflaton 
field. 
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According to a widely adopted picture [1], the perturbations to the Robertson- 
Walker cosmology arose from the quantum fluctuations in a slowly rolling 
scalar "inflaton" field during a period of inflation, then became classical as 
their wavelength was stretched beyond the horizon, and subsequently were 
imprinted on the decay products of the inflaton during a period of "reheat- 
ing." One of the attractions of this picture (and in particular the assumption 
that there is just one inflaton field) is that it has generally been thought to 
lead only to adiabatic perturbations, in agreement with current observations 
[2]. 

A recent preprint [3] has raised the question, whether it is possible for 
non-adiabatic cosmological perturbations to arise during reheating even af- 
ter inflation with a single inflaton field. This would be very important if 
true, for then observational limits on non-adiabatic fluctuations in the cos- 
mic microwave background would provide some constraints on the otherwise 
mysterious era of reheating, and indeed on the whole history of the universe 
between inflation and the present. 

However, there are very general grounds for expecting that single-fleld in- 
flation can only produce adiabatic fluctuations, whatever happens in reheat- 
ing or subsequently. It has been shown [4] that, whatever the constituents 
of the universe, the differential equations for cosmological perturbations in 
Newtonian gauge always have an adiabatic solution for wavelengths outside 
the horizon (that is, for physical wave numbers that are much less than the 
cosmological expansion rate). For this solution, there is a conserved quantity 
C [5], deflned by 

C = -* + <^p/3(p + p), (1) 

where 5 pis the perturbation to the total energy density in Newtonian gauge; 
bars denote unperturbed quantities; the perturbed metric is given by 

ds^ = -{I + 2^)dt^ + {1-2^) d^^ ■ (2) 

and as usual H = d/a. (Reference [4] dealt mostly with a quantity [6] 

n=-'^ + H5u. (3) 

where 6u is the perturbation to the total velocity potential, but outside 
the horizon TZ and C, are the same.) Also, for this adiabatic mode, the 
perturbations to the metric and the total energy density and pressure are 
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given outside the horizon by 



$ = ^ = ( -1 + 



a{t') dt' 



(4) 




5p^ 
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a{t') dt' , 
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This is a solution for any value of ti, so the difference of adiabatic solu- 
tions with different values of ti is also a solution, with C = ^''^^ $ = ^ oc 
H{t)/a{t), Sp/p oc l/a(t), etc. We will adjust ti so that the total pertur- 
bation takes the form (4), (5), in which case ti will be at some early time 
during inflation to avoid perturbations that are very large at early times. 
Furthermore, if the total energy-momentum tensor T^^ is given by a sum of 
tensors T^" for fluids labeled a (not necessarily individually conserved) then 
Eq. (5) holds for each of the individual perturbations 



and similarly for any other four-dimensional scalar, such as the inflaton field 



In general there may be other solutions, for which ( may or may not 
be constant and the perturbations are not given by Eqs. (4)-(7). (Mislead- 
ingly, these others are often called isocurvature solutions.) But if only one 
inflaton fleld makes a contribution to the energy-momentum tensor during 
inflation then the perturbations to this fleld and the metric in Newtonian 
gauge are governed by a third-order set of differential equations with a single 
constraint, so only two independent solutions contribute to cosmological per- 
turbations, and since we have found two explicit adiabatic solutions outside 
the horizon, these are the only solutions for the coupled system of inflaton 
and gravitational flelds from horizon exit to the end of inflation. During 
reheating and the subsequent evolution of the universe other fields and fiu- 
ids become important, but the adiabatic mode is a true solution outside the 
horizon whatever the contents of the universe, so with the universe in the 
adiabatic mode at the beginning of reheating it remains in this mode as long 
as the wavelength remains outside the horizon. 






C0(t) 



(7) 
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So how did non-adiabatic modes arise in reference [3]? The theorem of 
reference [4] assumes that the field equations hold at every moment, which 
requires that the perturbations are differentiable functions of time. As he 
suggested might be the case, the reason that Armendariz- Picon found a non- 
adiabatic mode during reheating in reference [3] was that, although it was 
assumed that initially during inflation only the inflaton flcid was present and 
there was no energy transfer to other fields, he supposed in this reference 
that the energy transfer rate rose discontinuously to a non-zero value at the 
beginning of reheating. (The model considered in reference [3] actually gave 
a constant value for but the perturbations had unequal values of 5p/p for 
the inflaton and its decay products, so as recognized by Armendariz-Picon, 
this perturbation was not adiabatic.) Of course, a discontinuous change in 
the energy transfer rate is unphysical. 

There is one weak point in the above general argument that non-adiabatic 
perturbations do not arise during reheating. It is the assumption that there 
was nothing but the inflaton and gravitation before reheating. Of course, 
in order for reheating to occur at all, there must be other fields or fiuids 
besides the inflaton, and these do not suddenly come into existence during 
reheating. The reactions that produce matter during reheating can't be com- 
pletely absent beforehand, so the remaining question is whether the transfer 
of energy from the inflaton to other fields during inflation before reheating 
excites these other modes. 

This question is answered by a further theorem, that if the matter en- 
ergy density during inflation is small, then even if the perturbations to the 
matter energy density are initially not at all adiabatic, the departures from 
adiabaticity would decay exponentially fast as soon as the energy transfer 
rate becomes appreciable. 

Here is the proof. The co-moving rate per proper volume of energy trans- 
fer from the inflaton field </? to "matter" fields (possibly including radiation) 
is in general some four- scalar function X of all these fields and perhaps their 
first and higher time derivatives: 

-uXm.^^^X{^,...), (8) 

where is the four-vector velocity of the total energy-momentum tensor, 
normalized so that u^u^ = — 1. To zeroth order in all perturbations to the 
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metric and other fields, this reads 

+ ^H{pm + Pm) , (9) 

where bars denote unperturbed values, taken to depend only on time. (We 
choose signs so that — To first order in perturbations, Eq. (8) gives 
in Newtonian gauge 

5pM + ^H{5pM + Spm) - 3(pm +Pm)^ ^5X + ^X , (10) 

the last term on the right arising from the perturbation to Uq. We have 
dropped terms involving spatial gradients, which become negligible outside 
the horizon. 

Now suppose that at some early time during inflation the density and 
pressure of matter are small. (This is plausible, because the energy density 
of fermions and gauge fields produced by quantum fiuctuations would be 
quadratic in the fluctuations.) Then the inflaton will provide the chief source 
of the gravitational fleld. As mentioned above, under these conditions the 
perturbations to the inflaton and gravitational potentials will be described 
by the adiabatic solution 

Sip = -ipx , $ = ^ = -c + i/x , (11) 

in which for convenience we have introduced the notation 

X{t)^^fa{t')dt'. (12) 

a{t) Jti 

Also, under these conditions the energy transfer rate will depend only on the 
inflaton and gravitational fields and their time derivatives, so 

SX = -XI . (13) 

(For instance, if X depends only on (p, then 

Eq. (13) also holds if X depends also on 0, p, etc., provided that for each 
pair of time derivatives there is a factor of g^^ to keep X a scalar.) Putting 
Eq. (13) together with Eq. (11), the right-hand side of Eq. (10) is 

d [ - ^ 

5X^^X^ -— XI . (14) 
at L J 
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and the difference between Eq. (10) and tlie time-derivative of Eq. (9) gives 



-Q^[hM + Pm'A = -^H{5pM + SpM)-HpM+PM)i-^-Q^ (pm+Pm)HX 

(15) 

Tflis can be simplified by noting tliat X satisfies tfie differential equation 

T+^(i/2:) = 0, (16) 

so 

d ._ ^ ^ 

[5pM + Pm^] = SH{5pM + SpM + Pm^ + Pm^) ■ (17) 

To continue, let us assume that the matter pressure is a function Pm{Pm) 
only of the matter energy density pm, as in the case of pure radiation or 
pure dust. (We are not assuming this for the combined system of matter and 
infiaton.) This is plausible, because the decay of a single real infiaton field 
would not generally produce any chemical potentials. Then pj^ — p'm{Pm)Pm 
and 8pM — p'm{Pm)8pm- Eq. (17) now reads 

d 

^ [SpM + PmI] = -^H{1 + cli){5pM + Pm'I) , (18) 

where c\,j = dpM/dpM is the squared sound speed. This shows that \5pM + 
Pm^] decreases monotonically and faster than a"^, but that is not good 
enough, because wc have to make sure that this is not just because \5pm\ 
and \Pm'^\ arc both decreasing. For this purpose, we use Eq. (9) again, and 
re- write Eq. (18) as 

^lnAr = -xii±iil, (19) 
ot Pm + Pm 

where is a dimensionless measure of the departure of the matter energy 
density perturbation from its adiabatic value 



N 



SpM + Pm^ 



pM +PM 



(20) 



Since energy is flowing from the infiaton field to matter, X is positive, and 
Eq. (19) shows that decreases monotonically. Early in inflation the matter 
perturbation may be nowhere near adiabatic, with N of the same order of 
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magnitude as the fractional density perturbation. As time passes during 
inflation, the transfer of energy from the inflaton to matter may make both 
5pM and pjy^ large, but the quantity N continues to decrease. Eventually, 
after the energy transfer rate X becomes large for a sufficiently long time 
during reheating, the matter density perturbation can no longer be ignored 
in calculating X and the gravitational field perturbations, but by that time 
will have decayed exponentially. With the density and pressure perturbations 
of the matter as well as the inflaton satisfying Eq. (6), we will still have 
5X — —XX and ^ — ^ — —X, and the above analysis will remain vahd. 

Note that asymptotically, after the transfer of energy from the inflaton to 
matter ceases, the quantity —pjyjX/{pM+PM) approaches — H/ H'^), so 

as long as pM < Pm, ~Pm^/ {Pm +Pm) is bounded below in absolute value by 
3|C|/4. Thus the reason that approaches zero is no^ that Spm/{pm + Pm) 
and —Pj^X/{pm + Pm) are both approaching zero. Rather, the ratio of these 
quantities, —Spm/Pm^i approaches unity. 

With the matter pressure a function only of the matter density, the matter 
pressure perturbation will then also satisfy Eq. (6). Thus the perturbations 
become adiabatic, in the sense of Eq. (6) and (7). Also, with the energy 
density and pressure perturbations of matter as well as of the inflaton all 
satisfying Eq. (6), the total energy density and pressure perturbations will 
obviously satisfy Eq. (5). It follows that the perturbations also satisfy the 
adiabatic condition that C, is constant, because in general 

The same analysis also obviously apphes if the inflaton energy goes into 
several species of matter, with a pressure for each species given by a function 
of the energy density in that species. 

In conclusion, even if the decay of the inflaton during inflation produces a 
small matter density whose perturbations are not at all adiabatic, the depar- 
ture from adiabaticity will decay rather than grow as inflation proceeds, and 
the departures of the perturbations from their adiabatic values will become 
exponentially small when the matter density becomes large during reheating. 

I was greatly helped in my thinking about this question by a correspon- 
dence some months ago with Alan Guth, and more recently by discussions 
with Christian Armendariz- Picon and Eiichiro Komatsu. This research was 
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